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Abstract 

We derive the asymptotic behaviour of the one point probability density for the 
inhomogeneous shock slopes in the turbulent regime, when a Gaussian fluctuating flux 
at origin derives the system. We also calculate the time dependence of the Xf beyond 
which there won't exist any velocity shocks as Xf = t 3//4 (log t) 1 / 4 . We argue that the 
stationary state of the problem would be equivalent with the long time limit of the 
diffusion equation with a random source at origin. 
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1 Introduction 



Burgers equation QTJ] initially was suggested as a model for describing the propagation of 
nonlinear waves in dissipative media |2[|. Burgers equation with a stochastic source has been 
suggested as a phenomenological equation covering the behaviour of a large universality class 
of nonequilibrium growth models ||. Under a nonlinear transformation it can also be trans- 
formed into a model for diffusion in Random landscapes which is itself related to a rich class 
of phenomena like, Directed Polymers in Random media [[|, wetting ||, Random Magnets 
H and Pinning of vortex lines in superconductors |7]. This equation is thoroughly studied 
as an approximate model describing the formation of large scale structures in Universe || . 
Burgers Equation with a noise term has also been studied as a toy model for fully-developed 
turbulence with using a variety of methods 0. The full description of this equation is still 
far from clear but the introduction of new methods has given a better understanding of the 
relationship between different realms of physics. The large time solutions of deterministic 
Burgers equation with random initial conditions has been investigated by the turbulence 



community for mimicking the behaviour of decaying turbulence ||10|| . An alternative way of 
introducing fluctuations is to consider random initial conditions. The physical relevance for 
instance comes from considering a variable injection of fluid into a semi-infinite slab of soil. 

Semiline solutions of Burgers equation with a given deterministic flux at origin has been 
derived |TTJ, which consistently reproduces the inhomogeneous shocks of Burgers equation 



at large spatial scales. By allowing random flux , we observe that shocks are present at large 
times but with random values and the velocities become ordered. 

2 The Semiline Solution 

Consider the Burgers equationon the semiline, with the initial and boundary condition: 

Ut + uu x = vu xx (1) 

u{x,0)=0 0<£<oo (2) 

^M + uu x (0,t)=g(t) (3) 
The above initial value problem is solved through a generalised Hopf-Cole transformation 

u(x,t) — —=—— — -j — T, — - — -, (4) 
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where: t 

c(t) = c(p)eJ° a W* f (5) 

In the limit where v tends to zero we can expand equation^ for regions where ^7= ^> 1: 

2 '° ^^SS) ea; ^-^4(fio " So' 9(t")dt")) 
u(x,t) = — - — , (6) 

^l + l^dt'g(t%^^)exp(-l(^-^g(tndtn) 
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Following Burgers 0, we note that the main contribution to the integral in (^j) come from 
points which maximize the exponent: 

h{x,tJ) = -^- T) + j\{f)df (7) 

After some manipulations, equation (|J) reduces to 

where ti(x,t) are the absolute maxima of h(x,t,t'). These are indeed the shock waves 
discovered by Burgers through an elegant geometric construction. 

Burgers geometric construction finds the maxima of h(x,t,t') by defining two functions: 

G(t') = /%(()«, (10) 

JO 

where if is a constant. For a given t, one looks for the maximum value of H, called H*(ti), 
where the hyperbola S(t') just touches the curve G(t'). For large values of x the hyperbola 
S(t') flattens, and it may touch G(t) at two points. When this happens we have a shock wave, 
thus the solutions are comprised of linear parts, rzir ■ In small regions of space, intervening 
the linear parts, the velocity profile is more complicated. The long time velocity profile is a 
series of shocks with reducing slopes. This can be understood as response of the system to 
impulses imported to it at origin, at consecutive times. 



3 One Point Probability Density 

Let us now analyse the same problem for the case of the random impulse at origin and 
assume that g(t) is a random variable with Gaussian distribution: 

(g(t)) = o (n) 

(g(t)g(t')} = 2a g 2 A(t-t') (12) 

Where A is any reasonable correlation which decays with growing argument. This means that 
G(t), as defined in equation (pi]) is a stochastic process and we assume that the variance of this 
process depends linearly on time so that it a nonhomogeneous process with some correlation 
time t overwhich its values are correlated and beyond that uncorrelated. Giving address to 
the geometrical construcion discussed in the previous section the problem reduces to finding 
the maxima of a nonhomogeneous process G(t). For this purpose we can adopt the methods 
previously used for the random initial conditions 0]. Consider the commulant probability 
of the event that the absolute maximum of h(x,t,t') occurs in the interval , t 2 ] ,that is, 

F(H, [t[,t' 2 ]) = Prob(-^—- + G{Q <H-Xe [t[,t' 2 ]) (13) 
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If the correlation time of G(t) is smaller than the mean difference of two adjacent ^ 's: 



1 



(t> +1 -t' i )>— 2 (14) 

then one can assume independence for regions inside the interval [ti, £2] arid outside it. Thus 
the probability of finding t t e [i' l5 1' 2 ] is : 

/oo 
F(fT, [0,t] - [t^M^, [0,t] - [4,4])^ (15) 
-00 



where P = dH 



dF 
dH' 

Indeed in the Poissonian approximation we have : 

FiH^tuh]) ^ e - N ^^ (16) 

where iV is the expected number of crossings between G(tf) and H + ^^tj in the interval 
[ti, t 2 ] and is given 

N(H,[t 1 ,t 2 ]) = (jm{^} + dt') (17) 

where <3> = ^_ t ,^ —G(f) and {x} + = x9(x). The average is taken over the joint distribution 
of g(t) and G(t), and there is the constraint that at the point of intersection G(t) has the 
positive slope. After some manipulations and in the limit when -J= <C t and H ^> a g t we 
get: 

2^ ,2 2 



N(H, [hM) = I ~ ^777^= e ( 18 ) 
Jti H\f2-K 



and consequently [0, t]) will become: 



Jo Sft' 
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F(#, [0,t]) = e Jo s (t ')v^ e (19) 

where S(t') = H + 4 ^_ t ,\ ■ For calculating the exponent we use from the saddle point 
approximation in the integral of the exponent. So in the large time limit we simply have the 
following saddle point solution for t* : 



x 2 



t -t*9i — + J— (20) 
AH \ AH y ' 

Before proceeding further we wish to find an estimate for Xf where the shock would disappear. 
Taking a closer look at the t* we find that when N(H, [0,t]) is maximum around t* the 
cummulant probability has it's most probable value and the contact point t* would be the 
most probable contact point. So a good estimate for Xf would be found by setting the 

t* = 0,from which we get 

x f VttH (21) 
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We believe that the main contribution for the integral (P^D would be around a character- 
istic value of if, say H* such that the cummulative probability F(H, [0, t}) gives its dominant 
contribution. On the other hand we find H* from 

N(H,[0,t]) S 1 (22) 

or 

a Q H' _ S(t') 2 

' =e 2tV « 2 = 1 (23) 



5(f) V2tt 

Intuitively this characteristic value gives the border beyond which there would be just one 
intersection between the process G(t') and the curve S(t'). After doing some algebra we find 
H* in the large time limit and in the domain i C as 



H*^a g Jt\ogt (24) 



By substituting the derived H* from (p4|) in expression fl2T|) we find the time dependence of 

*/ 2 = a// 2 (logt) 1/2 (25) 
When t \/fjf or i < ij the cummulative probability density will be: 

F(H,[0,t}) = e ^ ( * +a # (26) 

The asymptotic form of F(H, [0, t\) may be computed for large values of H = H* + f3z and 
the result comes out to be the well known Gumbel distribution \Y1 . 



F(H,[0,t})^e- e - z (27) 

where 

By substituting(|TB|) in ([T|) an taking an integration by part we will get 

P(t'e[t' 1 t' + At'}))= r N(H, [t t , h + At}) dF{H :Jt 0, t]) dH (29) 

JO (in 

We have calculated the probability distribution in (|29|) by approximating N(H, [£i, ii + At]) 
in fllBf ) for deviations from H* and substituting in fl2"9"|): 



P{t' e [f, t' + At'] \x, t) = 7 r(a/3 + 1, 1) (30) 
Where a and 7 are defined as 



>2o Q H' 1 (H * , * 2 , ^ 2 

: x_^ e ^(^+i(^) (31) 



5 



We have the following intuitive picture that as one measures the slope in a linear part of the 
velocity profile in a fixed spatial point the most probable value of is a decreasing function 
of x. By computing the numerical integration of the probability density one can easily see 
that for x's greater Xf which is given by (|2l| ) the total probability of occuring a contact is 
vanishing. By changing variable from t' to u(x, t) , one can easily find the probability density 
of u(x, t) for a fixed x and t as: 

P(u\x, t) = ^2a g ^-l)e-^^ {H 4 } r(^^ + 1, 1) (33) 

We have numerically integrated the average shock's slope and velocity for different values of 
x in a fixed observation time, and have plotted their character in different graphs (figures 
1,2). 



4 Two Point Probability Density 

In this section we generalise the above idea for calculating the two point statistics of the 
shocks. By using the geometrical picture of the contacts between the hyperbola and the 
random process Git'), every information about the two point statistics can be mapped to 
the different cathegories of geometrical events. If we are interested in the two point statistics 
in two different times ti, t 2 but at the same fixed point x, the correspondent event would be 
the contact of two parabolas with different vertical asymptotic axes at t\,t 2 , with the random 
process G(t') in two points say t[, t 2 without ever intersecting the other parts of it. Similarly 
the two point statistics in two different spatial points X\, x 2 ,and at a fixed observation time 
t, would be translated to definite geometrical events. Technically the same structure for 
calculating these quantities can be applied in principle. For example in calculating the two 
point probability density P(u\,u 2 ) on which u\ = u(xi,t) and u 2 = u(x 2 ,t), we use from 
the geometrical construction as in (fig 4). According to the figure, there are two situations 
distinguished which are in correspondence to existence and nonexistence of shocks between 
the two contact points t[ and t 2 . So we proceed to find 

P(t[ e U. A e L 2 ) = P p (t[,t' 2 ) + <*(<i - QPaiAA) (34) 

where P p corresponds to the peresence and P a corresponds to absence of shocks between two 
points t[ and t 2 . Existence of Dirac delta function in fl3~4] ) is related to the fact that in large 
spatial scales x G [x\ , x 2 ] and in the absence of discontinuties between t[ and t 2 the two 
hyperbolas should approximately contact near one of the local maxima of G(t), thus t[ = t 2 . 

By reffering to (fig 3) and adapting the techniques of previous section for calculating 
P p (t' 1 ,t / 2 ) we write: 

^ c -N(Hx,\t' +At,i\) d C -N{H U \H,,H+A^) d c -N(H 2 ,\t' A+At\) ^ 

dHi dH 2 
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on which t* is the intersection point of two hyperbolas, that is: 



w^h (36) 



It is clear that in the limit a g t 3> H* and for large times, the main contribution in is 
related to region near H* given in (E4j). Thus by defining Hi = H* + (3zi and H 2 = H* + /3z2 



the integration of (pop will give the two point statistics . The main point in calculating 
the integral is the existence of the same lavel H* over which the extremes of the proposed 
process G(t) have their major contribution. Calculations in this direction will be pursued in 
future. 



5 Diffusion Limit 

One of the interesting questions is related to stationary state of the system and its relation- 
ship to our picture. We believe that the stationary regime in our problem is governed by 
the diffusion equation stirred at origin. For showing the validity of our statement we should 
compare the average of the two terms in the denominator of the expression for u(x, t) given 
in (|4]). Before taking an integration by parts the denominator can be written as 

C (0)er/(-f=) + f dt> C{t ' ] e-^Fy (37) 

After taking the average over the distribution of g(t) one can easily find that the average of 
c{t) is 

(c(t)> = c(0)e^ (38) 

By substituting from (j38"|) into the ( j37|) and changing the variable of integration from i! to 
w we will end up with 



c(0)[er/(-=) + \ dww~^e-^ )w ~^ m ^\ (39) 



Because we are interested in steady state, we should find the large time limit of the denom- 
inator. The limit of the integral can be written in terms of the modified Bessel function 
Ki [4=] , so the denominator in large time limit can be written as 



C (0) e ^ (1_ ^ ) (40) 

where 

a H 2 

= ^ <«> 

It is obvious that in the limit where t — > +oo the denominator behaves as (^) which means 
that the integral in the denominator of ([|) would be negligible in average as if u(x, t) satisfies 
the diffusion equation with random source at origin. 
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6 Conclusion 



We have found the one point statistics of the Burgers equation with random boundary con- 
dition on a semiline. The random boundary condition in the problem can be translated as 
a random momentum input so the total momentum is not conserved. We have solved the 
problem by a geometrical construction which relates the original problem to the Extreme 
universality of the input noise in time(|T2|, |13[|). In the case of a non- homogeneous noise with 
a variance proportional with t we found the Gumbel universality class([|l^]). The overall pic- 
ture which we have found is that in a snapshot of the velocity profile there is a moving point 
Xdiff — y/vt* which gives a crossover region from the diffusive behaviour to shock structure 
such that with the passage of the time the diffusive region will grow and in the infinite time 
limit the diffusive regime would be dominant in the whole semiline. 

Recently it has been shown that there are some evidences which relates different Extreme 
universalities with different Replica symmetry breaking schemes (fl"3|, 14]) and specifically the 



relation between the Gumbel universality with to the one step Replica symmetry breaking 
is checked in some models. Checking another paradigm we think that it would be important 
to apply the Replica method on this problem too. Another problem worthy of investigation 
is the calculation of intermittency exponent structure function, and checking whether the 
nonhomogenity alters the intermittency exponent. Finally we think that the KPZ and di- 
rected polymer analogs of our problem which have been studied recently (]T5|) is relevant too 
because the known results in those realms are studied under the homogeneous noise and the 
presence of a nonhomogeneous noise might change the universality classes in those problems. 
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Figure Captions 

Figure 1- The behaviour of average shock slope in terms of x. a g — 1 and the observation 
time is t — 1000. The best fit for the curve is (A) = a+b ^, 3/2 where a = 182 and b = 0.05. 
Figure 2- The average velocity in terms of x.a g = 1 and the observation time is t — 1000. 
The best fit for the curve is In (u(x)) = a + b lj f where a = -0.1 and b = -16. 
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